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Abstract 

Core-collapse supernovae are among the most energetic cosmic cataclysms. They are prodigious emitters of neu- 
trinos and quite likely strong galactic sources of gravitational waves. Observation of both neutrinos and gravitational 
waves from the next galactic or near extragalactic core-collapse supernova will yield a wealth of information on the 
explosion mechanism, but also on the structure and angular momentum of the progenitor star, and on aspects of 
fundamental physics such as the equation of state of nuclear matter at high densities and low entropies. In this con- 
tribution to the proceedings of the Neutrino 2012 conference, we summarize recent progress made in the theoretical 
understanding and modeling of core-collapse supernovae. In this, our emphasis is on multi-dimensional processes 
involved in the explosion mechanism such as neutrino-driven convection and the standing accretion shock instability. 
As an example of how supernova neutrinos can be used to probe fundamental physics, we discuss how the rise time 
of the electron antineutrino flux observed in detectors can be used to probe the neutrino mass hierarchy. Finally, we 
lay out aspects of the neutrino and gravitational- wave signature of core-collapse supernovae and discuss the power of 
combined analysis of neutrino and gravitational wave data from the next galactic core-collapse supernova. 
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1. Introduction 

Core-collapse supemovae (CCSNe), primarily due to 
their nucleosynthetic impact on the chemical evolution 
of the universe, are of great importance in astrophysics. 
Their inner workings or, more precisely, the physical 
processes that convert the necessary fraction of gravita- 
tional energy of collapse into energy of the explosion, 
involve a broad range of physics, in particular neutrino, 
nuclear, and gravitational physics. The cores of CCSNe 
are thus cosmic laboratories for fundamental physics 
that can be probed by the observation and analysis of 
neutrinos and gravitational-wave signals from galactic 
or near-extragalactic CCSNe. This was first and im- 
pressively demonstrated by the observation of neutrinos 
from SN 1987A | V,T\, which confirmed the very basics 
of CCSN theory [3|: CCSNe are driven by the release 
of gravitational energy in the collapse of a massive star's 



core and the vast majority (~99%) of this energy is lib- 
erated in MeV neutrinos. 

CCSN theory and modeling have come a long way 
since the first pioneering numerical studies carried out 
in the late 1960s and early 1970s S IH EJ. Ad- 
vances in nuclear and neutrino physics, in combination 
with improvements in numerical modeling methodol- 
ogy and rapid advances in computer technology have 
since enabled detailed and self-consistent spherically- 
symmetric (ID), axisymmetric (2D), and the first (3D) 
models of CCSNe (see |[8l|3[ia for reviews). 

It is now clear that the prompt hydrodynamic mecha- 
nism proposed early on 1 1 1 1 fails robustly. The shock 
created at core bounce always stalls due to dissocia- 
tion of infalling heavy nuclei and neutrino losses from 
the region behind the shock. It must be revived by the 
explosion mechanism. The (delayed) neutrino mech- 
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anism lfT2l is now believed to drive most CCSN ex- 
plosions ||8]. It relies on the deposition behind the 
shock (resulting in net energy gain) of a small frac- 
tion (~10%) of the outgoing and neutrino lu- 
minosity within a few hundred milliseconds after core 
bounce. The neutrino mechanism has been shown to 
work robustly in spherical symmetry for the lowest- 
mass massive stars (e.g., fT3l). Driving explosions in 
more massive stars with the neutrino mechanism ap- 
pears to require an interplay of neutrino heating with 
multi-D fluid dynamics due to neutrino-driven convec- 
tion and the standing accretion shock instability (SASI) 
(e.g., ||ll|T5l[l6l[l7l[l8l[l9l|2a|2ll|22|23llMl and ref- 
erences therein). Convection/SASI leave characteristic 
imprints on the neutrino and gravitational-wave signals, 
which, in turn, can be used to observationally probe 
these instabilities l25l l26l l27l l28l l29l 1301 . 

While garden-variety CCSN explosions (of spectral 
type II, lb, Ic) may very well be ultimately explained 
by the neutrino mechanism, neutrino-driven explosions 
appear unfit to deliver hyper-energetic explosions (with 
energies of > 10 B; 1 [B]ethe - 10^' erg). Such hy- 
pernovae are predominantly of the rare Ic-bl subclass 
(bl stands for relativistically Doppler broadened spectral 
lines) of CCSNe that has also been associated with long- 
duration gamma-ray bursts (e.g., fT\\). Such hyperener- 
getic CCSN explosions may result from the magnetoro- 
tational mechanism. This mechanism relies on rapid 
progenitor star rotation and efficient amplification of 
protoneutron star magnetic fields to magnetar strength 
(~10'^G) after bounce, launching bipolar jet-like out- 
flows (e.g., fST, "331). The extreme rotation required in 
this mechanism may leave characteristic imprints on the 
neutrino and gravitational wave signal Il34ll35l . 

An interesting potential twist to the standard CCSN 
story developed with the (re-)discovery of collective 
neutrino oscillation induced by neutrino-neutrino co- 
herent forward scattering (see, e.g., |36| for a review). 
Supernova neutrinos, as reviewed by Balantekin in this 
volume, may oscillate in vacuum, due to the MSW ef- 
fect ||37l , and, in regions of high neutrino density, due 
to collective oscillations. In the inverted neutrino mass 
hierarchy, collective oscillations are predicted to lead to 
a nearly complete swap of Ve/v^ spectra with the spectra 
of heavy lepton (x - fi, r) neutrinos vjv^. If this swap 
occurs below the region in which neutrino heating dom- 
inates over neutrino cooling, this would potentially lead 
to a boost of neutrino energy deposition, since the emis- 
sion spectra of Vx/v^ are significantly harder than those 
of their /v^, counterparts and the charged-current ab- 
sorption cross sections scale with the square of the neu- 
trino energy (e^). That this could indeed have an im- 



portant impact on neutrino-driven explosions was first 
shown by Suwa et al. [38]. However, recent detailed os- 
cillation calculations on the basis of realistic supernova 
background matter distributions and neutrino radiation 
fields found that collective oscillations are likely to be 
suppressed in the early post-core-bounce evolution and 
thus cannot enhance neutrino heating ll39l l40l 1411 (but 
see also (42]). 

In this contribution to the proceedings of the Neu- 
trino 2012 conferences, we outline, in Section |2] some 
recent results from multi-D simulations of CCSNe. In 
Section [3j we discuss how the early postbounce CCSN 
neutrino signal can be used to constrain the neutrino 
mass hierarchy. Finally, in Section [4] we give examples 
of what may be learned from combined observations of 
neutrinos and gravitational waves from the next nearby 
CCSN. 

2. Results from Recent Neutrino-Driven Multi- 
Dimensional Core-Collapse Supernova Models 

The first set of detailed 2D CCSN simulations carried 
out in the mid 1990s (e.g.. El |44l |43) demonstrated 
that the negative entropy gradient created by neutrino 
heating leads to vigorous convection in the gain region 
behind the shock (where neutrino heating dominates 
over cooling). The added degree of freedom in 2D liber- 
ates the infalling gas to move not just in the radial direc- 
tion (as in ID), but also laterally. As a consequence, the 
dwell time in the gain region of a given accreting mass 
element is increased, leading to greater total neutrino 
heating and high-entropy turbulent flow that pushes the 
shock to larger radii and aids the explosion. 

The standing accretion shock instability (SASI; e.g., 
Il46l l47l l48l ) is another hydrodynamic instability of 
relevance in the postbounce evolution of CCSNe. In 
2D simulations, it drives large-scale periodic low-order 
(i = 1 is the fastest growing mode) deformations of 
the shock front, which, in the nonlinear phase of the in- 
stabiUty, lead to larger average shock radii, secondary 
shocks, and also increase the dwell time of material in 
the gain region (e.g.. Ell HI [IB [HI). In 3D, non- 
axisymmetric (\m\ > 0) modes will develop. If rotation 
is present, m + I spiral modes may appear (e.g., |49l). 
while in the nonrotating case all m of a given { are de- 
generate and should grow to comparable amplitudes, ef- 
fectively spreading the typical 2D { amplitude over its 
m-modes in 3D |50, 22, 21. 23 J. 

There has been much recent discussion in the litera- 
ture about the respective roles of neutrino-driven con- 
vection and SASI and about which of the two instabil- 
ities is the primary instability (i.e., the one that grows 
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Figure 1: 3D volume rendering of tiie specific entropy i at ~ 150 ms 
after core bounce in a 3D general-relativistic simulation caiTied out 
by Ott et al. 1231 . The scale of the frame is 700 km on a side. 
The colormap is chosen such that cyan correesponds to a moder- 
ate « ~ 4 kB baryon" ' . Yellow and red indicate regions with i ~ 
16kB baryon"' and s ~ 20 kB baryon"', respectively. 

first and is dominant) |[l6l[l7l[ll|22l|2T]|2Ql|2l. While 
this is not the place to review this discussion in detail, it 
is worth pointing out a few key aspects: 

(/) Which instability grows first depends on (a) the 
early postbounce rate of advection through the region 
of convective instability (and, hence, on the accretion 
rate set by progenitor star structure) and (b) on the size 
of the seed perturbations present in the flow from which 
convection may grow p48', "471 and thus on the level of 
asphericity in the last stages of nuclear burning (e.g., 
iBTI ). If the initial size of a perturbation is too small to 
develop into a bouyant plume before it is advected out 
of the region of convective instability, then convection 
cannot grow and SASI will dominate (at least initially) 

Eiiini. 

(ii) When the SASI approaches its nonlinear phase, 
it leads to secondary shocks and seeds convection 
(e.g., pTl), which, in turn, may be limiting further SASI 
growth (e.g., |52|). 

(///) The behavior of the SASI appears fundamen- 
tally different in 3D and 2D. In 2D, the dominant £ - 
1 mode leads to strong periodic north-south sloshing 
motions of the shock that persist even when neutrino- 
driven convection becomes strong. In 3D, ^ = 1 is 
still the fastest growing mode. However, its power is 
equipartitioned into its m - {0, ± 1 ) modes, the invidi- 



3 

ual ({, m) amplitudes stay at much smaller values than 
in 2D and no large-scale single sloshing mode develops 
Il50l[l8l|29l|22l|2l]|23]. 

(iv) In 3D, if neutrino-driven convection gets a chance 
to grow, it appears to robustly become the dominant in- 
stability and quench the SASI, in particular in simula- 
tions that have strong neutrino heating and show ex- 
plosions II22I l2n I23I . In the latter, non-oscillatory, 
high-amplitude { = l,m = {0, ± 1 ) deformations ap- 
pear that correspond to large bubbles of high-entropy 
that push out the shock front in a particular direction 
Il22l|2l]|23|53l. An example of this is shown in Fig.[l] 
which depicts a volume rendering of the specific en- 
tropy at ~ 150 ms after core bounce in a simulation of 
Ott ef a/.ll35l. 

Nature is 3D. While highly detailed 2D multi-group 
radiation-hydrodynamics simulations are now available 
lfT6HT7l . a full understanding of CCSN hydrodynamics 
must await equally detailed 3D simulations, which have 
not been performed, but rapid progress towards them is 
being made (e.g., l54ll55lE3l). 

3. Supernova Neutrinos as Probes 
of the Neutrino Mass Hierarchy 

Recently, Serpico et al. If56l proposed that the neu- 
trino mass hierarchy could be determined from the early 
postbounce neutrino signal of the next galactic CCSN 
for which complicated collective oscillation effects may 
be ignored ||39ll40| . The original flavor (at emission) of 
Vg arriving at a detector on Earth depends on the neu- 
trino mass hierarchy. For the normal hierarchy (NH), 
the MSW efi'ect ll37l in the outer regions of the star leads 
to an incomplete conversion of Vg to heavy-lepton an- 
tineutrinos (vj) and vice versa. For this hierarchy, in the 
adiabatic approximation for neutrino oscillations, the Vg 
signal at Earth is a linear combination of the original Vg 
signal (oc cos (612)^ ~ 0.71) and the original signal 
(oc sin(0i2)^ ~ 0.29). However, for the inverted hierar- 
chy (IH) the Vg signal at Earth is completely composed 
of the original signal 157J . So, for the two hierarchies, 
we have: 

L^^^ = COS(012)%„SN + sin(0i2)%„SN , (1) 
^9^,e - ^v„SN • (2) 

The distinguishing characteristic explored in ll56l . 
and further in this proceedings contribution, is the dif- 
ference in the rise times of the Vg signal and the Vi sig- 
nal. The production of Vg is initially suppressed due 
to the large electron chemical potential before bounce 
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Figure 2: Normalized (to the luminosity at 100 ms after bounce) elec- 
tron antineutrino (Vj) luminosities at Earth for 32 progenitor star mod- 
els and assuming both the inverted (family of curves labeled via Vj) 
and normal (family of curves labeled O.VlVf +0.29;i"« v) hierarchy. The 
colors distinguish progenitor models, colors ranging from red — > yel- 
low — ► green — ► blue denote increasing compactness, Eq. (|3j 1581 . 

and in the very early postbounce. It is not until the de- 
generacy has been partially lifted in the outer regions 
of the protoneutron star that significant luminosity 
emerges from the supernova. Hence, Vg have a slow rise 
time. Heavy-lepton neutrino production does not un- 
dergo this suppression and the Vj luminosity rises much 
more quickly. In the following, we explore whether a 
high-event-rate detector like IceCube (e.g., Il59l ) will be 
able to distinguish between the two expected neutrino 
signals given by eqs. ffli and (j2]i. In particular, we ex- 
tend the analysis of 1561 and explore whether there ex- 
ist degeneracies in the observed signal with progenitor 
star model and equation of state (EOS). We use results 
of O'Connor & Ott ||58][ ^ who simulate the preexplo- 
sion phase in 32 solar-metallicity progenitor stars with 
masses ranging from 12 to 120 Mq from Ref. |60|. Two 
equations of state are used, the LS220 EOS from 1611 
and the HShen EOS from f62]. In Fig. [S] we show 
the normalized (to its value at 100 ms after bounce ) 
luminosities expected at Earth assuming both the nor- 
mal (slowly rising curves) and inverted hierarchies (fast 
rising curves). We show only the luminosities for the 
LS220 EOS. The fast rise time of the inverted hierarchy 
is clearly visible in the data. This figure is very similar 
to Fig. 4 of 1 56 1, but shows a larger range of progenitors, 
which we parameterize by their compactness parameter, 

/;(Mba,y = M)/1000kmlr=Wc ' 



Figure 3: Ratio of the number of expected events detected in the first 
20 ms to the number of expected events detected in the first 100 ms 
versus cumulative events at 100 ms. We show this ratio for each pro- 
genitor star model, EOS, and most important for our discussion, the 
assumed neutrino mass hierarchy. The fast rise time of the inverted hi- 
erarchy signal gives a large ratio for all models, independent of EOS 
and progenitor star details. Note that neither axis begins at zero. Er- 
ror bars are purely statistical and include the contribution from the 
IceCube background rate of ~2200 events per 1 .6384 ms bin. 



where RiMh^^ - M) is the radial coordinate that en- 
closes a baryonic mass of M at the time of core bounce. 
Here, we use M - 1.75 Mq, since this is a fiducial mass 
scale in the postbounce preexplosion phase [58 1. For 
clarity, the color scale in Fig. |2] ranging from red 
yellow — > green blue, denotes increasing progenitor 
star compactness ll58l . 



' Available at jhttp : //www . stellar collapse ■ org/Mlprog 



We predict interaction rates in IceCube by using the 
following formula \35 \, Rf^ ~ 2.6 x lO'^s^'L® jj/r^gj^^^, 
where L? is the v^, luminosity at Earth in units of 
10^' erg s"' and riokpc is the distances in units of 10 kpc, 
the fiducial galactic distance. We quantitatively dis- 
criminate between the two signals at Earth by taking 
the fluence ratio of the number of events in the first 
20 ms to the number of events in the first 100 ms. While 
we have not performed an analysis to find the optimal 
times to evaluate the fluence ratio, 20ms/100ms is a 
reasonable choice, as one can see in Fig.|2] It is appar- 
ent that by 20 ms the Vi luminosity has begun to level 
off. If we picked a time smaller than 20 ms, we would 
start to be dominated by statistical error for the num- 
ber of observed events. Times in the numerator longer 
than 20 ms lead to fluence ratios that are increasingly 
converging towards one. The inverted hierarchy, which 
gives a faster rising signal in at Earth, will have a 
larger fluence ratio than the normal hierarchy. We show 
this ratio in Fig. |3] for each progenitor model, EOS, and 
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hierarchy. For clarity, we choose to plot the fluence ra- 
tio as a function of total number of event^ detected at 
100 ms. The results of our extended model set confirm 
the results of obtained by ||56l . The normal and inverted 
hierarchies should be distinguishable from the rise time 
of the early neutrino signal from the next galactic CCSN 
in an IceCube-like detector. In addition, over the range 
of progenitors that we explore and the two quite distinct 
EOS we employ, the signals predicted by the two hier- 
archies remain clearly distinct, strengthening the ability 
to use supernovae neutrino rise times as a probe of the 
neutrino mass hierarchy. 

4. The Next Galactic Core-Collapse Supernova: 
Gravitational Waves and Neutrinos Probing the 
Supernova Engine 

The next galactic CCSN is guaranteed to come. 
While the predicted rates for galactic CCSNe are in 
the range of a disheartening one per ~30 - 100 years 
ll63l l64l . we prefer to look at this from a more posi- 
tive perspective: The next galactic CCSN has already 
exploded. Using the predicted range of rates, one 
can, using a simple toy model of the galaxy assum- 
ing uniform star formation in the disk, estimate that 
-1000 (SNrate)/(3.33x 10-2 SNe yr-^) CCSNe have al- 
ready exploded and their photons, neutrinos, and gravi- 
tational waves (GWs) are on their way to Earth. The ar- 
rival of these messengers from the next galactic CCSN 
may be imminent. Due to intervening dust and gas, a 
significant fraction, perhaps 30 - 50%, of these CCSNe 
will not be discovered in the optical, but will rather be 
seen as very luminous infrared sources and/or via their 
X- and y-ray emission. Photons observed from CCSNe 
will initially carry information on the outer layers of the 
progenitor star, and only at late times allow us to probe 
central regions. Neutrinos and GW, on the other hand, 
will provide plentiful information on the thermodynam- 
ics and dynamics of the supernova engine. 

As shown by Lund et al. |28 | (also seen in ll65l[34l 
IT5I ). the SASI leads to characteristic time variations in 
the neutrino flux at ~50 - 100 Hz that are observable by 
large Cherenkov detectors (see |66| and other contribu- 
tions in this volume for a discussion of supernova neu- 
trino detectors) for a galactic CCSN. If neutrino-driven 



^We note that unlike the ratio of events at 20 ms to 100 ms, the to- 
tal number of events detected depends on the distance to the CCSN, 
which may not be determined (see However, given that the hierar- 
chies are distinguishable for all models regardless of the total number 
of detected events (i.e. progenitor and EOS), this will not limit their 
distinguishability. 



convection dominates, the preexplosion time variations 
of the neutrino flux may be expected to be characterized 
by smaller amplitude and higher frequency. Hence, the 
neutrino signal of the next galactic CCSN may obser- 
vationally constrain the prominence of neutrino-driven 
convection and SASI. 

GWs probe accelerated aspherical (of lowest-order 
quadrupole) mass-energy motions. Detailed reviews of 
the CCSN GW signature can be found in |l25]|67l. In 
the case of a nonrotating or slowly rotating CCSN, GW 
emission will occur via (/) prompt postbounce convec- 
tion, which may grow on the negative entropy gradi- 
ent left behind by the stalling shock, (//) convection 
inside the protoneutron star, (///) convection/SASI be- 
tween the protoneutron star and the stalled shock, (;v) 
asymmetric neutrino emission, and (v) aspherical out- 
flows. The upcoming second-generation of km-scale 
GW interferometers (Advanced LIGO |68J, GEO-HF 
Il69l . Advanced Virgo CO), and KAGRA lEIl) and 
third-generation instruments such as the Einstein Tele- 
scope 172] will be most sensitive to GWs in the fre- 
quency band ~30 - 1000 Hz, and will most likely detect 
GW emission from galactic and near extragalactic CC- 
SNe (within, perhaps, ~100kpc). Emission processes 
{i-iii) lead to broadband emission in ~50 - 1000 Hz, 
while iiv) and (v) are characterized by low-frequency 
(< 10 Hz) GWs and are thus less relevant for Earth- 
based detectors. 

Studies investigating GW emission from convec- 
tion/SASI (e.g., 122, 26i|23l|29l|23i^) have shown that 
the strongest component of GW emission comes from 
the deceleration of undershooting downflow plumes at 
the edge of the protoneutron star. As the region of de- 
celeration is changing in the early postbounce phase, 
the time-frequency evolution of the GW signal observed 
from the next galactic CCSN will provide information 
on the early evolution of the structure, composition and 
thermodynamics of the lower postshock region, set, pri- 
marily, by the nuclear EOS and the postbounce accre- 
tion rate ll27l [30l . This information encoded in GWs 
could be combined with information from neutrinos, 
which also probe the structure of the outer protoneutron 
star 1581. 

In the case of a rapidly spinning CCSN progenitor 
that may explode via the magnetorotational mechanism, 
core collapse proceeds approximately axisymmetrically 
11741 175]| . At bounce, the extreme acceleration experi- 
enced by the rotationally-deformed {t - 2, m = 0) in- 
ner core leads to a strong burst of GWs, followed by 
quasi-periodic ring-down oscillations which continue 
for ~10 - 15 ms. At later postbounce times, rotational 
nonaxisymmetric (m + 0) dynamics may develop due 
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to a shear instability and lead to longer-term quasi- 
periodic GW emission at typically twice the rotation 
frequency of the protoneutron star Il74ll75l . The charac- 
teristics of the GW signal due to dynamically-relevant 
rotation are rather distinct from those of the GW signal 
due to convection/SASI. If it is true that rapid rotation 
leads to magnetorotational explosions, while slowly ro- 
tating/nonrotating progenitor stars explode via the neu- 
trino mechanism, the GW signal alone may be sufficient 
to determine the explosion mechanism via Bayesian 
model selection \76\. 

The neutrino signal seen from a rapidly rotating 
CCSN wiU depend on the inclination of the spin axis 
with respect to the observer (e.g., (34] and references 
therein): The rotational flattening of the protoneutron 
star leads to higher neutrino fluxes with harder spectra 
along the poles (i.e. the rotation axis) and to lower fluxes 
with softer spectra in equatorial regions. The hardness 
of the neutrino spectrum is also sensitive to the nuclear 
EOS (e.g., [38" and the inclination of the CCSN 
will not be known, unless it goes along with a gamma- 
ray burst, indicating a rotation axes aligned with the line 
of sight. Additional information carried in the GW sig- 
nal can be used to break the degeneracy between ro- 
tation and EOS and constrain the unknown inclination 
angle. For example, if rotation is rapid, the protoneu- 
tron star ring-down oscillations seen in the GW signal 
go along with time variations in the Ve and Vi fluxes at 
the same frequency that may be detected out to 0(1 kpc) 
distances |35 1. If the CCSN spin axis is pointed at Earth, 
the {( - 2,m - 0) GWs emitted by the ring-down oscil- 
lations will be suppressed, but the time-variations in the 
neutrino fluxes will be seen, as will be the subsequent 
GW emission by nonaxisymmetric dynamics, which is 
of ({ - 2,m = 2) character and peaks along the rota- 
tion axis f75\. Combining these different emissions will 
allow observers to constrain the inclination of the mag- 
nitude and orientation of the CCSN spin, in addition to 
helping to break the degeneracy between rotation and 
nuclear EOS. 

The potential of combined analysis of GW and neu- 
trino signals from the next galactic CCSN is tremen- 
dous. Much theoretical and modehng work is still 
needed to flesh out the details and parameter depen- 
dence (progenitor structure, rotation, EOS etc.) of both, 
GW and neutrino signals. This will require further im- 
provements in modeling technology and efliciency to 
enable parameter studies in 3D. Progress is also needed 
in the development of parameter estimation approaches 
that will allow us to make statistical statements based 
on the combination of information carried in GWs and 
neutrinos. 
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